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ABSTRACT 

We present the first non-LTE time-dependent radiative-transfer simulations of supernovae 
(SNe) II-Plateau (II-P) covering both the photospheric and nebular phases, from ~10 to 
>1000d after the explosion, and based on 1.2 B piston-driven ejecta produced from a 15 M Q 
and a 25 M Q non-rotating solar-metallicity star. The radial expansion of the gradually cooling 
photosphere gives rise to a near-constant luminosity up to ;>100d after the explosion. The 
photosphere remains in the outer 0.5 M Q of the ejecta for up to ^50 d after the explosion. As 
the photosphere reaches the edge of the helium core, the SN luminosity drops by an amount 
mitigated by the progenitor radius and the 56 Ni mass. Synthetic light-curves exhibit a bell- 
shape morphology, evolving faster for more compact progenitors, and with an earlier peak and 
narrower width in bluer filters. UV and <7-band fluxes are very sensitive to line-blanketing, the 
metallicity, and the adopted model atoms. During the recombination epoch synthetic spectra 
are dominated by H I and metal lines, and are largely insensitive to the differing H/He/C/N/O 
composition of our two progenitor stars. In contrast, synthetic nebular-phase spectra reveal 
a broader/stronger O I doublet line in the higher-mass progenitor model, reflecting the larger 
masses of oxygen and nickel that are ejected. Our simulations overestimate the typical lu- 
minosity and the visual rise time of standard SNe II-P, most likely a consequence of our 
progenitor stars being too big and/or too hydrogen rich. Comparison of our simulations with 
photospheric-phase observations of SN1999em of the same color are satisfactory. Our neglect 
of non-thermal excitation/ionisation leads to a fast disappearance of continuum radiation and 
Balmer-line emission at the end of the plateau phase. With the exception of H I lines, our 
nebular spectra show a striking similarity to contemporaneous observations of SN1999em. 

Key words: radiative transfer - stars: atmospheres - stars: supernovae - stars: supernovae: 
individual: SN 1999em 



1 INTRODUCTION gen: IWooslev & Weayedri995t) . as well as efficient dust factories 

JSugerman et all2006h . 

Supe rnovae (SNe) II-Plateau (II-P) are the mos t represented SN 

type dCappellaro etafl 1 19971 ; ISmartt et alj l2009h . The theoretical ° win g to their large luminosity, SNe II-P offer a means to con- 

prediction th at they stem from the explosion of Red Su pergiant strain distances in the Universe through a variety of methods known 
(RSG) stars iColgate & Wfutell 19661 : lFalk& Arnettlll977l) is sup- as the Expanding Photosphere Metho d (iKirshner & Kwarj|l974l: 
ported by direct identification of some SNe II-P's with their pro- lEastman et al.ll 19961; iDessart & Hilliej|2005all2006h, the Spe ctral- 
genitor on pre-explosion images. The progeni tors indicate main- fitting Expanding Atmosphere Method Mitchell et alj|2002l). and 
sequence masses in the range 8-17 M (see ISmartJ [20091 for a the Standard Candle Method faamuv & Pintol2002l). Studying SNe 
review). After their death, SN II-P produce a stellar -mass com- n " p is thus of relevance for a wide range of topics in modern astro- 
pact object, either a neutron star or a black hole dHeger et all physics. 

2003). They are a key driver for the chemical evolution of galaxies Inferences on SNe II-P ejecta, as well as their progenitor stars, 

through the core-embedded material they eject (in particular oxy- can be made using radiation hydrodynamics and radiative trans- 
fer. For multi-band light-curve calculations, radiation hydrodynam- 
ics assuming the gas in Local Thermodynamic Equilibrium (LTE), 
frequency-mean or frequency-dependent opacities (including lines 
* E-mail: Luc.Dessart@oamp.fr or not) and time-dependent radiation transport are used with some 
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success to constrain the explosion energy, the mass of the ejecta, 



in the explosion (e.g., 


Falk & Arnettll 19771: llitvinova & Nadezhinl 


19851; Blinnikov et al] 


200d: iBaklanov et al.ll2005l; Utrobinll 20071: 


Tominaga et al.l 12009; 


Dessart et al. 2010b. aS). Calibrations drawn 



Table 1. Properties of the pre-SN models sl5 and s25 iWoosley & Heger 
2007). The age is the time elapsed between the main sequence and core 
collapse. 



from such stu dies have been used to infer SNe II-P properties 
dHamuy||2003l) . 

However, an attractive and more physically-consistent alter- 
native is to use radiation hydrodynamics up to a few days after 
shock-breakout, when acceleration terms are sizeable, and then 
switch to a more sophis ticated radiative-tran s fer treatment on a 
non-accelerating ejecta ([Eastman et alj 1 19941 ; iKasen & Wooslevl 
l2009l : iDessart & Hillienl201Ch . This then permits a more detailed 
treatment of line-blanketing, departures from LTE, and the depen- 
dence of the radiation field on frequency, angle, and time. Re- 
cently u _w£_sketched_aii_a£P_roach that treats all these important as- 
pects toessart&HillieJl2010h and applied it to the early evolution 
of SN1987A. The method retains the sophistication of non-LTE 
steady- state simulations focusin g on the decoupling layers of SN 
ejecta dDessart & Hillied l2005bl ; iBaron et alj 120071 ; I5essart et al.l 
2008), but also treats important time-dependent terms appearing in 
the energy, statist ical-equilibrium, and radiative-transfer equations 
(see, for example, lUtrobin & Chugail2005l ; IDessart & Hillied 2008 , 
for a partial treatment of these terms). 

This initial-value problem is constrained at the start by em- 
ploying a hydrodynamical input model of the explosion, itself gen- 
erated from a physically-consistent pre-SN model. Confrontation 
of synthetic spectra and light curves with observations can thus 
be connected to the known ejecta properties, yielding constraints 
on the explosion and the progenitor properties. Because we ex- 
plicitly treat all processes controlling level populations, the ef- 
fects of line emission, absorption, and overlap are accurately com- 
puted, and thus line profiles can be used to extract information 
on the ejecta. Accurate spectral modelling is critical for the de- 
termination of abundances, excitation/ionisation sources, kinemat- 
ics, and ultimately the explosion and progenitor properties. By con- 
trast, most SN studies treat bound-bound transitions approximately. 
These studies allow inferences based on colour and light-curve evo- 
lution but little reliable information is extracted from lines. 

In this paper, bu ilding upon our recent analysis of the Type 
II-peculiar SN1987A dDessart & Hillierll2010h . we present the first 
non-LTE time-de pendent radiative-transfer simulations of SNe II- 
P usin g CMFGEN dffillier & Milledl 19981 ; IDessart & Hillieill2005bl 
2008), covering both the photospheric and nebular phases. Be- 
cause of the huge computational expense of this endeavour, we fo- 
cus on only two ejecta produced from the explosion of stars with 
main-sequence masses of 15 Mq and 25 Mq, and characterised by 
an explosion energy of 1.2 B (Woosley, private communication; 
IWoosley & Hegeij|2007l) . We use these simulations to describe fun- 
damental properties of the SN gas and light, with special emphasis 
on spectra which have not previously been computed with this level 
of physical consistency in this context. Due to the limited set of in- 
puts, we make no attempt to reproduce specific observations, but 
rather we seek to make a qualitative assessment of the general ade- 
quacy of our models. 

In the next section, we present the pre-SN models that were 
used to generate the hydrodynamical input ejecta from which 
our radiative-transfer simulations start. The adopted model atoms, 
atomic data, and the importance of non-LTE, are discussed in §[3] 
In §[4] we discuss the evolution of the gas properties throughout the 
photospheric and nebular phases, including the photospheric prop- 
erties and the ejecta ionisation. We then turn to our theoretical re- 
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suits for the radiation field. In §[5] we first describe the bolometric- 
light evolution and multi-band light-curves, before moving on to 
synthetic spectra, the sources of blanketing, and line identifications 
in §|6] We also document the spectral differences predicted in our 
simulations for the low and the high-mass progenitor models. In 
§ [7] we compare our results to photometric and spectroscopic ob- 
servations of SNe II-P in general, and of SN1999em in particular, 
to identify the successes and failures of our simulations. In § [8] 
we present our conclusions and plans for a forthcoming parameter 
study of SNe II-P. 



2 INITIAL HYDRODYNAMICAL MODEL AND SETUP 

The approach we f ollow is similar to that described in 
IDessart & Hillied d2010l) . but we now concentrate on the SN ejecta 
resulting from the explosion of RSG stars. Starting from a hy- 
drodynamical model that describes the depth variation of ele- 
ment mass fraction, density, temperature, radius and velocity, we 
model the non-LTE tim e-dependent radiative t r ansfer of SNe II- 
P ejec t a with CMFG EN dHillier & Milled 1 19981 ; IDessart & Hillieil 
2005b, I2008L |2010|) . We assume the SN is spherically symmet- 
ric and smooth (i.e., un-clumped), and free from external distur- 
bances, such as interaction with the circumstellar environment or 
irradiation from the newl y-born neutron star. The h ydrodynamical 
models were computed by Woosley & H eger ( 2007) using the code 
KEPLER dWeaver et alJI 19781) by first evolving non-rotating, solar- 
metallicity, 15 Mq and 25 Mq stars from the main sequence to the 
formation and collapse of their degenerate iron core, and then ex- 
ploding these by means of a piston. 

The inpu t model s!5e!2 (s25e!2) was generated from model 
sl5A (s25) of IWoosley & Heged d2007l) and corresponds to a SN 
ejecta with a 1.2 B energy and post-explosion time of 8xl0 5 s 
(2x 10 6 s; see Table Q]for details). For the same ejecta kinetic en- 
ergy, the model s25el2 yields a higher mass of 56 Ni (0.163 com- 
pared to 0.0866 Mq) because it possesses a larger helium-core 
mass with more mass at larger density, smaller radii and higher 
temperatures. All models were artificially mixed by sweeping a 
"box" four times through the ejecta and averaging the composi- 
tion within that box. The box used in model sl5el2 (s25el2) had 
a width of 0.427 Mq (0.819Mq), corresponding in each case to 
a tenth of the progenitor helium-core mass. This moderate mixing 
softens the composition gradients but preserves the strong ejecta 
stratification, hydrogen being absent in the inner few solar masses 
of the ejecta (below 1500 km s" 1 ). The 56 Ni mass fraction is size- 
able only below 1000-2000 km s~ x but owing to mixing, shows a 
non-zero mass fraction even at the outer edge of the ejecta. 

In the (standard) KEPLER models sl5el2 and s25el2, the 
species included are H, 3 He, 4 He, 12 C, 14 N, 16 0, 20 O, 24 Mg, 2S Si, 
32 S, 36 Ar, 40 Ca, 44 Ti, 48 Cr, 52 Fe, 54 Fe, 56 Ni, 56 Fe, and "Fe". In 
this nomenclature, some species approximate certain nuclei. For 
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Table 2. Ejecta cumulative yields (upper table) and surface mass fraction (lower table) for the most abundant species in models sl5el2 and s25el2. Metal 
abundances not shown here are taken at the solar-metallicity value. The total ejecta mass is 10.93 Mq (13.94Mq) for model sl5el2 (s25el2). Model sl5el2iso 
has the same ejecta characteristics as model sl5el2 and differs only in composition of under-abundant species, e.g., Sc or Ba. Numbers in parentheses 
correspond to powers of ten. 



Model 














Ejecta Cumulative Yields 














A/ H 

[Mq] 


M He 

[Mq] 


M C 

[Mq] 


M N 

[Mq] 


M 

[Mq] 


M Ne 

[Mq] 


M Mg Af Si M s 

[Mq] [Mq] [Mq] 


Mat M Ca 

[Mq] [Mq] 


M Ti 

[Mq] 


M Cr 

[Mq] 


M Pe 
[Mq] 


Af56 Ni 
[Mq] 


sl5el2 
s25el2 


5.49(0) 
4.01(0) 


4.00(0) 
4.77(0) 


1.59C-1) 3.22(-2) 
3.53(-l) 4.48(-2) 


8.16(-1) 
3.32(0) 


1.52(-1) 
5.30(-l) 


5.41(-2) 8.10(-2) 3.17(-2) 
2.24(-l) 3.16(-1) 1.49(-1) 


6.43(-3) 5.52(-3) 
2.58(-2) 1.85(-2) 


1.26(-4) 
1.54(-4) 


1.29(-3) 
2.48(-3) 


1.60(-2) 
2.04(-2) 


8.66(-2) 
1.63(-1) 


Model 












Ejecta Surface Mass Fractions 
















X C 




Xo 




^Mg -^Si X§ 


X\ r Xca 






^Fe 


X&6 Ni 


sl5el2 6.50(-l) 
s25el2 5.28(-l) 


3.35(-l) 
4.56(-l) 


1.42(-3) 
9.17(-4) 


3.11(-3) 
5.54(-3) 


5.41(-3) 
3.89(-3) 


1.31(-3) 
1.37(-3) 


7.91(-4) 8.29(-4) 4.23(-4) 
8.03(-4) 8.33(-4) 4.25(-4) 


1.13(-4) 7.38(-5) 
1.13(-4) 7.40(-5) 


3.83(-6) 3.43(-5) 
3.83(-6) 3.43(-5) 
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Figure 1. Left: Element mass fractions as a function of Lagrangian mass for the hydrodynamical input model sl5el2. We show hydrogen (black), helium 
(pink), oxygen (turquoise), carbon (blue), nitrogen (green), and 5e Ni (red). We also overplot the velocity as a dotted black line, with the corresponding 
axis given on the right. We include an axis at the top showing the distribution of the mass-integrated specific kinetic energy (the integral bounds being the 
Lagrangian mass M r and Aftot). The bulk of the ejecta kinetic energy is contained within what used to be the loosely-bound hydrogen-rich region of the 
progenitor envelope. In contrast, the more gravitationally bound helium core is ultimately endowed with only a few percent of the total ejecta kinetic energy. 
Right: Same as left, but now for the s25el2 model. The helium core is endowed with a larger fraction of the total ejecta kinetic energy than seen in the sl5el2 
model, although the largest share is still contained in the hydrogen-rich region of the ejecta. For both models, the total ejecta kinetic energy is 1.2B. 



example, 4 He approximates nuclei with 1 < A < 6, where A is 
the atomic mass, and 48 Cr stands for nuclei with 47 < A < 52 
[see Woosley, Heger, & Weaver (2002) for details]. In such inputs, 
numerous species are not described individually (e.g., Na, although 
Na I lines are unambiguously see in SNe II-P spectra and require 
the knowledge of the sodium abundance for a proper modelling). 
Some, like Ba, are altogether ignored. We give a summary of the 
chemical composition for the dominant species in Table [2] with 
both the cumulative ejecta yields (upper table) and the correspond- 
ing mass fractions at the progenitor surface (lower table). We also 
illustrate the distribution of the main elements through the ejecta 
in Fig. Q] Surface abundances of both models show evidence for 
CNO processing, with the higher-mass model exhibiting larger en- 
hancements of both He and N, and a greater depletion of C and O. 
The presence of CNO processed material in RSG a tmospheres is 
corro borated by spectroscopic analyses of RSG stars dLancon et al.l 
120071) . and understood to stem, in the absence of rotation, from con- 



vective dredge-up. Stellar rotation is expected to aid the appearance 
of CNO processed materia l at the stellar surface ( Mevnet & Maederl 
2000: iMevnet et alj |2006). In the outer regions of the ejecta, mass 
fractions for elements heavier than oxygen are identical for both 
models and correspond to the solar metallicity (apart from the nu- 
cleosynthesised 56 Ni). 

The bulk of the kinetic energy is contained within the 
hydrogen-rich regions of the ejecta, while the region correspond- 
ing to the helium core contains <10% of the total kinetic energy 
in model s25el2, and an even lower fraction in the model sl5el2 
(see top axis in Fig.QJ. The higher mass progenitor (which suffered 
enhanced mass loss) has a lower mass hydrogen-rich envelope and 
a larger helium/oxygen ejected mass (Fig. [TJ. As we discuss in this 
paper in relation to the degeneracy of photospheric-phase SNe II-P 
spectra and light curves, the enhanced mass of ejected oxygen may 
be most easily identified in SNe II-P observations taken at nebular 
times. Because it pertains to the systematic increase in helium-core 
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Figure 2. Top row: Comparison between the initial KEPLER hydrodynamical input sl5el2iso (black) and the corresponding ejecta structure we start with for 
the radiative-transfer simulations (red). The differences between the two arise from our need of strict-homologous expansion (r/v = constant) and of an 
outer optically-thin region where photons are freely-streaming (which requires we stitch this outer region, which is lacking, to the KEPLER model). In the 
process, we artificially speed up regions with v < 2000 km s _1 by <50% and smooth the density structure of the SN ejecta. We indicate in the middle panel 
the density exponent N p = —d log p/d log v for the inner, intermediate, and outer ejecta regions. As discussed in Section|4] the photosphere resides in the 
region of steep density decline with N p ~ 20 — 25 and 7.5 < Vph o t/(1000 km s~ 1 )< 9 for up to three weeks after explosion in model sl5el2iso. Bottom 
row: Same as top row, but now for model s25el2. The speed-up of the inner regions is smaller, approximately 10%, because we start with a model that is 
closer to homologous expansion. 



mas s with main-sequence , it can help constrain the progenitor iden- 
tity jPessart et all2010al) . 

In parallel to the simulation based on model sl5el2, we 
also perform a CMFGEN simulation based on the detailed KE- 
PLER model sl5el2iso, which has exactly the same structure in 
density/radius/ velocity/temperature, but was evolved with a more 
complete nuclear network. The KEPLER model sl5el2iso includes 
every element of the periodic table (as well as all its associated iso- 
topes) up to Astatine (symbol At; Z =85). This detailed input per- 
mits the specific modelling of any species (and associated ions). In 
practice, apart from unstable isotopes, we use the sum of the indi- 
vidual isotopic mass fractions to determine the abundance of the as- 
sociated species. Compared to models sl5el2 and s25el2, our sim- 
ulation of model sl5el2iso treats in addition the elements Ne, Na 
(no longer assumed to have an abundance at the solar-metallicity 
value and constant with depth), Al, Ar, Ti, Sc, Cr, and Ba, adopt- 
ing the abundance distributions computed by KEPLER. With this 
choice, we include all species that have a mass fraction larger than 
~10 at any depth in the sl5el2iso ejecta, as well as any species 
with identified features in SNe II-P spectra (e.g., Sc and Ba). When 
discussing general ejecta and radiation properties, we employ our 
simulations based on the sl5el2 and s25el2 models, but when 
comparing to observations we employ our simulations based on 
model sl5el2iso as it describes more completely the ejecta com- 
position and the effects of line blanketing. 

Before starting our radiative-transfer calculations, we made a 
number of adjustments to the hydrodynamical inputs. In both mod- 
els, the innermost regions are not in homologous expansion (despite 
the ejecta age of 9.26 and 23.15 d for models sl5el2 and s25el2, 
respectively) so we enforce r/v = constant (adopting a midpoint 
in the ejecta) for compatibility with the radiative-transfer algorithm 
which assumes dv/dr = v/r. This artificially enhances the small 
velocity of the inner regions by ^50% in model s!5el2 and by 



~10% in model s25el2. To ensure the radiation is free streaming at 
all wavelengths at the outer boundary, and to correct for the poorly- 
resolved outer region in the input model, we trim and stitch an outer 
region to the input model (characterized by a power-law density 
with an exponent of 15), extending it to larger radii/velocities. The 
photosphere is always located inside this outer region, and is thus 
not directly affected by this manipulation (see Section . At the 
new outer boundary, densities are on the order of 10 -20 gem -3 
and velocities on the order of 20000 km s . To illustrate these 
adjustments, we show in Fig. [2] the ejecta properties provided by 
KEPLER and those we start with in our CMFGEN simulations. 

One issue of concern is the starting conditions. In our ap- 
proach we fix the temperature to that of the hydrodynamical sim- 
ulation and solve the relativistic transfer equation and statistical 
equilibrium equations neglecting time dependent effects. For the 
models presented here this is likely to be a reasonable approxi- 
mation. Steady-state CMFGEN models showed that time-dependent 
effects for SN 1999em did not b ecome noticeable in spectra un- 
til day 16 jPessart & H illier 2006). On the other hand our starting 
times would be unreasonable for SN 198 7 A, where time dependent 
effects are already imp ortant at day 4 dUtrobin & Chugail 120051 : 
iDessart & HilTie3l2O10h . The reasons for starting these exploratory 
simulations late were two fold — the need to run many additional 
models at earlier time steps was considered unnecessarily wasteful, 
and the assumption of a homologous expansion at earlier times is 
increasingly problematic. 

The post-explosion time corresponding to our "adjusted" in- 
puts has been modified by 5-8% (because of the enforcement of 
homology) so the ejecta age is 8.5 1 d (22.01 d) for our initial simu- 
lation based on model sl5el2 (s25el2). To step in time, we adopt a 
fixed logarithmic increment, At/t ~ 0.1 and continue the simula- 
tions just past 1000 days (we also run simulations using a fixed 2-d 
time increment to cover the rapid fading at the end of the plateau 
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Time since Breakout [d] Time since Breakout [d] 

Figure 3. Evolution of the photospheric radius iiis (black; Ri$ = JJphot/10 15 cm), photospheric velocity V4 (blue; V4 = V p h o t/10 4 km s — 1 ), photo- 
spheric temperature T4 (red; T4 = T p ^ ot /W 4 K), and the mass above the photosphere AMi (green; A Mi = AMphot/10 Mq) until 150 d after explosion, 
and shown for simulations based on model sl5el2 (left) and s25el2 (right). The dots correspond to the actual times at which the radiation-transfer simulations 
were performed. The evolution at nebular times is excluded since no photosphere exists at such epochs. 
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sl5el2 (II-P) 
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1 10 
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100 



Figure 4. Evolution of the photospheric velocity obtained with our simu- 
lations based on SNe II-P models sl5el2 and s25el2 (blue and red lines; 
this work) and the SN II-pec model lm87Ael2 of SN1987A (black line; 
iDessart & Hilliej2010l) . 



phase). For the simulation based on model sl5el2 (s25el2), the se- 
quence covers until 1040 d (1037 d) in 52 (43) steps and took a bout 
3 months to complete for each. As in lDessart & Hillietl d2010l) the 
gas temperature is held fixed, and set to that provided by the hydro- 
dynamical input model, for the first model in the sequence. For sub- 
sequent models, the temperature is solved for everywhere. Hence, 
only results from the second time step and onwards are shown or 
discussed here (corresponding to ages 9.3 and 24.0 d for simula- 
tions based on models sl5el2 and s25el2). We impose a floor value 
to the gas temperature to avoid problems with the evaluation of ex- 
ponentials involving high ionisation stages. Since the ejecta cool as 
time proceeds, it eventually becomes unnecessary to include high- 
ionization stages. As we discard them, starting from the highest 
ones, we can lower the floor temperature and not face the floating- 
point-errors we would otherwise get in the evaluation of exponen- 
tials. In practice, we use a floor temperature of 2000-4000 K dur- 
ing the plateau phase, lowered to 1000-2000 K during the nebular 
phase. A lower floor temperature, or none at all, would be more 



suitable. Indeed, our tests suggest that this is in part the reason for 
the overestimated flux in the UV and U B bands. 

The spatial grid used for the radiative-transfer solution is 
set to have (approximately) a constant increment in optical-depth 
(computed using the Rosseland-mean opacity), with typically 5 
points per optical-depth decade, and a maximum of 100 points. 
This is much less than the original 750-900 data points of the 
hydrodynamical-input grid. However, our grid is much finer for 
the regions where radiation decouples from matter, the transition 
from Eddington factors of 1/3 to 1 occurs over a few tens of points, 
while in the hydro input this transition takes place over merely 
a few points. Furthermore, as the ejecta recombine, our grid ad- 
justs to track the recombination fronts (across which the optical 
depth varies steeply with radius). In SNe II, three main fronts form, 
associated with the recombination to Hell, He I, and Hi. How- 
ever, these front s never steepen as fier cely as in our simulations 
of SN 1987 A jDessart & Hilliedl2010h o r as in steady-state SNe 
II-P simulations (e.g.. |Pessart et alj|2008l) . Convergence is not in- 
hibited by this feature and is indeed quite steady in the simulations 
we describe here. 

For the treatment of radioactive decay, we assume only one 
decay chain corresponding to 56 Ni — > 56 Co — > 56 Fe. In a forth- 
coming paper, we will present a Monte Carlo code that solves for 
the transport of 7-ray photons from the site of emission out through 
the ejecta and potentially towards escape. We find that in SNe II, the 
deep-seated location of nucleosynthesised 56 Ni and the large over- 
lying mass prevent the escape of 7-ray photons up to >500 d. Non- 
local energy deposition starts however to operate at a few hundred 
days. As the ejecta become thin to optical photons near the end of 
the photospheric phase, high-energy (i.e., non-thermal) electrons, 
Compton-scattered by 7-rays, will affect the excitation and ionisa- 
tion level of the gas. 

A proper and accurate treatment of radioactive decay re- 
quires first the determination of the energy depositi on versus 
depth dAmbwani & SutherlandH"l988l; ISwartz et alj|l995l) . and sec- 
ond the heating/exci tation/ionisation effects of t h e high-energy 



nergy 

electrons produced jAxelrodl Il980l; iLucvl Il99ll; ISwartz! Il99lt 
i Kozma & Franssonlll992k ISwartz et al J 1 1 9931 : iKozma & Franssonl 
Il998allbl) . Although aware of these important issues, we assume in 
this first study of SNe II-P that all decay energy is deposited locally 
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as heat. We anticipate this will overestimate the SN visual bright- 
ness at nebular times, and underestimate the excitation/ionisation 
level of the gas. In our simulations, this neglect seems to affect 
primarily hydrogen, causing B aimer lines and the associated con- 
tinuum to vanish abruptly beyond ~130d after explosion, contrary 
to observations (see §[7}. 



3 MODEL ATOMS, ATOMIC DATA, AND DEPARTURES 
FROM LTE 

The model atoms used in our sim ulations are similar to those 
adopted in iDessart & Hillier! J20ld) . The choice is a compromise 
between physical accuracy, completeness, and computational bur- 
den. Our simulations require about 6GB of RAM and one time step 
takes about 2-3 days on a single CPU. In CMFGEN, we can in prin- 
ciple treat all species up to Z — 28 (Nickel), as well as Barium 
(Z=56). However, to limit the computational burden, our early CM- 
FGEN simulations based on models sl5el2/s25el2 neglected some 
species and/or ions. In Table[3] we document these details, includ- 
ing for each ion species the number of super and full levels, the 
number of transitions, and the uppermost level included^] To cap- 
ture the higher level of detail of the ejecta composition for model 
sl5el2iso and compute more accurate spectra for comparison with 
observations, we additionally include in our treatment Ne II, Al II, 
Al III, Ar III, Sc I, Sc II, Sc III, Cr II, Cr III, Fe I, Ni II, Ni III, 
Ni IV, Ba I, and Ba II (these additional ions appear with a super- 
script a). As we discuss below, our simulations based on sl5el2iso 
show the importance of treating Na I, Sc II, Cr II, and Fe I, as well 
as considering the corresponding elemental abundance stratifica- 
tion with depth. 

The model atoms described in Table[3]are used to start the time 
sequences. As time proceeds, the ejecta recombine and cool so that 
the highest ionisation species can progressively be discarded. Dur- 
ing the nebular phase, our model atom typically contains species 
that are at most three-times ionised. Concomitantly, the low den- 
sity and weaker radiation field at late times favours the appear- 
ance of forbidden lines (associated with Ca II and O I in particu- 
lar). We therefore adapt the super-level assignment to capture this 
change in conditions as ions get less and less excited. In practice, 
past 100-200 days, we typically split the lower 10 super levels for 
all species (splitting the lower 20 super levels does not alter the 
results visibly and supports this choice). The set of atomic pro- 
cesses, m ost of the atomic data, and their use, have been sum- 
marised in |Pessart& Hillier! i feOlOl) and discussed in more detail in 
iHillier & Millerl dl998t) ; they are thus not repeated here. For Ba II, 
no t included in earlier models, we used the atomic data discussed 
bv lMazzali et al.l d 19921) (kindly supplied by K. Butler), while fo r 
Sell and ScIII we used oscillator strengths from iKuruczl {20090 
and crude approximations for the photoionisation cross-sections. 
Energy levels, when available, were taken from NIST0 



1 Level designations are generally according to the National Institute of 
Standards and Technology (NIST), except we use "z" to refer to states 
where several high "1" states have been merged (say I = 4 to n — 1), 
and "w" states refer to a level where all "1" states for a given multiplicity 
have been merged. 

2 Data is available online at http://kurucz.harvard.edu 

3 Data is available online at http://physics.nist.gov/asd3 



Table 3. Summary of the complete model atom used in our radiative- 
transfer calculations for model sl5el2iso (our simulations for models 
sl5el2 and s25el2 were done without the species iden tified with super- 
script a). The source of the atomic datasets is given in Dessart & HillieJ 
120101) . Nf refers to the number of full levels, N s to the number of super 
levels, and Ntrans to the corresponding number of bound-bound transi- 
tions. The last column refers to the upper level for each ion treated. The 
total number of full levels treated is 8273, which corresponds to 173965 
bound-bound transitions. 



Species 


N f 


N s 


Ntrans 


Upper Level 


Hi 


30 


20 


435 


n ^ 30 


He 1 


51 


40 


374 


n ^ 11 


Hell 


30 


13 


435 


n ^ 30 


Cl 


26 


14 


120 


n ^ 2s2p 3 3 Pg 


Cn 


26 


14 


87 


n ^ 2s 2 4d 2 D 5/2 


Cm 


112 


62 


891 


n ^ 2s8f 1 F° 


Civ 


64 


59 


1446 


n ^ 30 


Ni 


104 


44 


855 


n ^ 5f 2 F° 


No 


41 


23 


144 


n s£ 2p 3 d 1 P 1 


Ol 


51 


19 


214 


n ^ 2s 2 2p 3 ( 4 S)4f 3 F 3 


On 


111 


30 


1157 


n ^2s 2 2p 2 ( 3 P)4d 2 D 5/2 


Om 


86 


50 


646 


n ^ 2p4f L D 


OIV 


72 


53 


835 


n ^ 2p 2 ( 3 P)3p 2 P 


OV 


78 


41 


523 


n ^2s5f *F° 


Nell a 


242 


42 


5794 


n ^ 2s 2 2p 4 !D4d 2 S 


Nal 


71 


22 


1614 


n SC 30w 2 W 


Mgll 


65 


22 


1452 


n ^ 30w 2 W 


AlII a 


44 


26 


171 


n ^ 3s5d X D 2 


AlIIII a 


45 


17 


362 


n ^ lOz 2 Z 


Sill 


59 


31 


354 


n s:3s 2 ( 1 S)7g 2 G 7/2 


Sim 


61 


33 


310 


n ^ 3s5g G4 


Si IV 


48 


37 


405 


n SC lOf 2 F° 


Sll 


324 


56 


8208 


n 3s3p 3 ( 5 S°)4p 6 P 


S in 


98 


48 


837 


n ^3s3p 2 ( 2 D)3d 3 P 


SIV 


67 


27 


396 


n Os3p( 3 P°)4p 2 D 5/2 


ArIIII a 


346 


32 


6898 


n sg 3s 2 3p 3 ( 2 D°)8s ^ 


Call 


77 


21 


1736 


n s; 3p 6 30w 2 W 


Scl a 


72 


26 


734 


n s£ 3d4s( 3 D)5s 2 D 5/2 


Scll a 


85 


38 


979 


n s£ 3p 6 3d4f 1 V 1 


Sclll" 1 


45 


33 


235 


n sS7h 2 H 11/2 


Till 


152 


37 


3134 


n s£ 3d 2 ( 3 F)5p 4 D 7/2 


Tim 


206 


33 


4735 


n s: 3d6f 3 Hg 


CrII a 


196 


28 


3629 


n s; 3d 4 ( 3 G)4px 4 G 11/2 


CrIII a 


145 


30 


2359 


n 3d 3 ( 2 D2)4p 3 Dg 
n 3d 6 ( 5 D)4s4p x 5 Fg 


Fe'I a 


136 


44 


1900 


Fe II 


1 15 


50 


1437 


n $ 3d u ( i Gl)4s d z G 7 / 2 


Fem 


477 


61 


6496 


n sg 3d 5 ( 4 F)5s 5 Fi 


FelV 


294 


5 1 


8068 


n sg 3d 4 ( 5 D)4d 4 G 5/ / 2 


Fe V 


191 


47 


3977 


n ^ 3d 3 ( 4 F)4d 5 F3 


Fevi 


433 


44 


14103 


n ^ 3p5( 2 P°)3d 4 ( 1 S) 2 P 3/2 


Fe VII 


153 


29 


1753 


n ^ 3p5( 2 P°)3d 3 ( 2 Dl) ^1 


Coll 


144 


34 


2088 


n ^ 3d 6 ( 5 D)4s4p 7 D° 


Com 


361 


37 


10937 


n ^ 3d 6 ( 5 D)5p 4 P 3/2 


Co IV 


314 


37 


8684 


n s£ 3d 5 ( 2 P)4p 3 Pi 


Cov 


387 


32 


13605 


n s£ 3d 4 ( 3 F)4d 2 H 9/2 


Co VI 


323 


23 


9608 


n s£ 3d 3 ( 2 D)4d 1 S 


Co VII 


319 


31 


9096 


n ^ 3p5( 2 P°)3d 4 ( 3 F) 2 D 3/2 


Nill a 


93 


19 


842 


n s; 3d 7 ( 4 F)4s4p 6 D° /2 


NilII a 


67 


15 


379 


n ^ 3d 7 ( 4 F)4p 3 D° 


NilV a 


200 


36 


4085 


ns£3d 6 ( 3 D)4p 2 D 5/2 


NiV 


183 


46 


3065 


n s£ 3d 5 ( 2 D3)4p 3 F§ 


Ni VI 


314 


37 


9569 


n s£ 3d 4 ( 5 D)4d 4 F 9/2 


Ni VII 


308 


37 


9225 


n s£ 3d 3 ( 2 D)4d 3 P 2 


Bal a 


31 


17 


30 


n s£ 5d6p ^1 


Ball a 


100 


38 


2514 


n ^ 3s30w 2 W 
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AM phot [M s ] AM phot [M ] 

0.01 0.10 2.27 6.16 9.72 10.67 0.05 0.26 4.48 11.63 13.47 




Time since Breakout [d] Time since Breakout [d] 

Figure 5. Evolution of the composition at the photosphere for a few important species in our simulations based on models sl5el2 (left) and s25el2 (right). 
The dots give the epoch for each simulation in the time sequence, and AMphot m Fig.[3]is the mass contained between the photosphere and the outer-model 
boundary. Notice the coincidence of the end of the plateau phase (shown further below in Fig. |5J with the recession of the photosphere into the helium-rich 
(hydrogen-poor) layers of the ejecta. 



3.1 Departure Coefficients 

In Fig. [8] we show the departure coefficients as a function of the 
Rosseland-mean optical depth (ra osse iand) for H I, O I, and Fe II at 
two different epochs. These departure coefficients are defined as the 
ratio of the non-LTE population to the LTE population, where the 
LTE population has been computed using the electron temperature, 
and the actual values of the ion and electron densities. As expected, 
LTE is recovered at large optical depth (rR oaae i an d> 10). Below 
7"Rosseiand=10, significant departures from LTE begin to occur, and 
these are well developed around the photospheric layers. At small 
optical depths, the departure from LTE can be orders of magnitude. 
In general the behavior of the departure coefficients is complex, and 
highly dependent on the level under consideration. 

The very significant departures of the populations from their 
LTE values is not surprising. In general electron scattering, and line 
scattering dominate the opacity, and these do not strongly couple 
the radiation field with the electron temperature. Further, the elec- 
tron density is low at the photosphere (~ 10 10 cm -3 ), and hence 
collisions are generally ineffective at causing a significant coupling 
of the majority of the level populations and ionisation state of the 
gas, and cannot overcome their stronger coupling with the radiation 
field. 

Above the photosphere the electron temperature is strongly 
decoupled from the radiation field temperature, and is instead de- 
termined by an equilibrium between many different processes, in- 
cluding adiabatic cooling, collisional cooling and heating, and heat- 
ing by photoionization. Indeed, at the outer boundary of sl5el2iso 
at, e.g., ~50d, the temperature has reached the floor value of 
3000 K, while the photospheric temperature is somewhat lower 
than 6000 K. If we assume that the ground-state population is de- 
termined by the radiation temperature, we would expect a depar- 
ture coefficient of order exp(157, 850/T rad - 157,850/T o ) or 
~ 3.0 x 10~ 12 assuming T e = 3000 K, T rad = 6000 K and a 
fixed electron density; hence very small departure coefficient for 
the lower levels are to be expected. Of course the situation is more 
complex than assumed above — the radiation field is a diluted 
black-body, and in the outer regions the H ionization is not in equi- 
librium because the recombination time-scale is longer than the ex- 
pansion time scale. 

For Hi, the departure coefficients are close to unity for high 



Time Since Breakout [d] 
10 267 525 782 1040 




2000 4000 6000 8000 10000 
Velocity [km s -1 ] 



Time Since Breakout [d] 
10 267 525 782 1040 




2000 4000 6000 8000 10000 
Velocity [km s -1 ] 

Figure 6. Top: Evolution of the ejecta gas temperature from 10 to 1040 d 
after explosion in the simulations based on model sl5el2. Dots give the 
position of the photosphere, which by definition exists as long as the base 
electron-scattering optical depth is greater than 2/3. A color coding is used 
to differentiate the epochs. Bottom: Same as top, but now showing the elec- 
tron density. 
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principal quantum number, as expected. For Fe II, all the levels 
shown have very small departure coefficients at low TR OSSe i an( j. 
This is not surprising — although we have shown 115 levels (with 
a smaller number of super levels), the principal quantum number of 
the valence electron for levels illustrated is less than 5. For Fe II the 
lower metastable levels are coupled by collisions with the ground 
state, and hence tend to be in LTE with respect to the ground state 
(i.e., they have the same departure coefficient as the ground state). 
This can be seen in Fig. [8] — many levels have similar departure 
coefficients around TR OS seiand=l. A similar statement applies to 
metastable levels of other elements. 

In Monte-Carlo simulations of SN spectra it is common prac- 
tice to use approximate formula to determine the ionization struc- 
ture, and to use an LTE assumption to dete rmine populations (e.g., 
Kasenetal . 2006l; lKasen & Woo slev 2009). As the photon scatter- 
ing is treated explicitly, the calculations are neither LTE, or rigor- 
ously non-LTE. As apparent from the above, and indeed from our 
need to explicitly allow for the time-dependence when solving for 
the hydrogen ionization structure, non-LTE populations need to be 
allowed for when modeling Type II SN. Experience indicates that 
as we include more processes, more lines, and more levels, the de- 
partures from LTE will probably become smaller, but the departures 
will still be highly significant. 



4 EJECTA EVOLUTION 

Over the 10-1000 d covered by each simulation, the ejectum ex- 
pands by two orders of magnitude and changes from being rather 
dense and optically-thick to being very tenuous and optically-thin. 
Over the same time period, the ejectum cools from a few 10 5 K 
at depth, to temperatures everywhere less than 3000 K. These dra- 
matic changes cause an equivalently dramatic change in the radi- 
ated light, which we describe in §[5]and §|6] 

During the photospheric phase, the electron- scattering radial 
optical depth is, by definition, greater than a few. Thus during this 
phase there always exists a "photosphere", which we formally de- 
fine as the location where that optical depth is 2/3. The behaviour of 
the photosphere determines to a large extent the SN light curve, and 
we illustrate the evolution of its properties in Fig.[3] Starting ;>10d 
past explosion, our simulations miss the phase of rapid evolution of 
photos pheric properties which immediately follows shock break- 
out dPessart et alj|2008l; iGezari et alJbood : iTominaga etafll2009l ; 
Rabinak & Wax manll2010l) . Indeed, at the start of our simulations, 
the photosphere has already expanded by about a factor of ten, its 
temperature dropped to <J10000K, and its velocity has been re- 
duced to <8500km s~ , Subsequently, for each simulation, and 
over the first ~80d, we find a slow decrease of the photospheric 
temperature and velocity, but a systematic increase in radius. 

Importantly, the photosphere resides within the outer 0.5 M© 
of the ejecta for about ~50d after the explosion (Fig.[5]l- Over that 
initial expansion, the ejecta has sufficient thermal energy to remain 
ionised while expanding, so that the photosphere remains close to 
the progenitor surface layers. As recombination sets in during the 
plateau phase, the photosphere eventually recedes in mass, but this 
recession is initially compensated by expansion so that the change 
in photospheric radius is small — in the present models the photo- 
spheric radius only begins to shrink after ~ 80 d. 

The photospheric composition does not change for the first 
~70 d in these two SNe II-P models an d is identical to that of the 
progenitor surface. The same finding bv lEastman et al.l d 19941) sug- 
gests that this is a generic property of SNe II-P. It implies that the 



Table 4. Summary of photospheric-phase properties for the non-LTE time- 
dependent simulations based on model sl5el2. For each epoch, we give 
the total ejecta electron-scattering optical depth, the photospheric radius (in 
10 15 cm), the photospheric velocity (in km s _1 ), the photospheric temper- 
ature (in K), the amount of mass above the photosphere (AAf p h t, given in 
Mq), as well as the emergent bolometric luminosity (in 10 s Lq). Results 
for the time-sequence based on sl5el2iso are the same as for those based 
on si 5e 12, and hence are not shown here. 



sl5el2 



Age 




^?phot 


Vphot 


7~phot 


AM phot 


Lbol 






Liu cmj 


[Km s | 


I'M 


UVlQ J 


rin8i ~\ 
LIU L0j 


9.30 


14034.20 


0.69 


8545 


10661 


0.008 


7.94 


10.20 


11053.89 


0.74 


8406 


10063 


0.010 


7.64 


11.20 


9134.01 


0.80 


8294 


9527 


0.011 


7.54 


12.30 


7567.96 


0.87 


8183 


9057 


0.014 


7.50 


13.50 


6259.40 


0.94 


8083 


8642 


0.016 


7.44 


14.80 


5192.15 


1.02 


7988 


8317 


0.019 


7.37 


16.30 


4252.72 


1.11 


7890 


8054 


0.023 


7.24 


17.90 


3490.37 


1.21 


7800 


7859 


0.027 


7.11 


19.70 


2823.76 


1.31 


7701 


7678 


0.032 


6.96 


21.70 


2287.51 


1.42 


7595 


7477 


0.038 


6.78 


23.90 


1863.47 


1.54 


7480 


7208 


0.046 


6.65 


26.30 


1519.30 


1.67 


7337 


6878 


0.057 


6.56 


28.90 


1233.58 


1.79 


7164 


6564 


0.075 


6.49 


31.80 


993.68 


1.91 


6957 


6305 


0.101 


6.42 


35.00 


796.48 


2.03 


6715 


6087 


0.141 


6.35 


38.50 


630.28 


2.14 


6438 


5893 


0.200 


6.30 


42.40 


495.07 


2.25 


6135 


5735 


0.290 


6.29 


46.60 


388.62 


2.32 


5760 


5478 


0.455 


6.19 


51.30 


304.74 


2.38 


5378 


5403 


0.733 


6.19 


56.50 


240.12 


2.44 


4999 


5355 


1.118 


6.15 


62.10 


190.81 


2.49 


4634 


5306 


1.619 


6.03 


68.30 


149.01 


2.52 


4277 


5255 


2.274 


5.85 


75.10 


113.56 


2.54 


3918 


5211 


3.097 


5.61 


82.60 


83.02 


2.54 


3556 


5174 


4.071 


5.32 


90.00 


58.83 


2.50 


3214 


5162 


5.069 


5.02 


99.00 


37.12 


2.43 


2835 


5185 


6.159 


4.52 


108.00 


21.80 


2.28 


2441 


5276 


7.144 


3.98 


118.00 


9.53 


1.97 


1932 


5415 


8.069 


3.31 


130.00 


2.19 


1.44 


1286 


5775 


9.723 


2.18 


143.00 


0.75 


0.86 


698 


5550 


10.628 


1.11 



spectral evolution of SNe II-P during the first half of the plateau 
phase does not stem from composition changes, but instead from 
changes in ionisation. The same ap plies to the early evolution of 
SN 1987A dPessart & Hilliej2010l) . 

The near uniform photospheric composition with time during 
the photospheric phase, as well as the modest composition con- 
trast between the h yd rogen envelopes of RSG models (see, e.g., 
IWooslev et al 1120021 or lWooslev & Heger 2007), highlights the dif- 
ficulty of constraining the progenitor identity with photospheric- 
phase spectra. It also partially explains the uniformity of SNe II- 
P spectra, which tend to show lines from the same ions, with the 
same strength and shap e s, and generally o nly differing in width 
dPessart & HillieJ 120061 ; IPessart et alj|2008l) . This suggests that 
differences between SNe II-P are difficult to extract after the onset 
of recombination in the photospheric phase. They may instead be 
constrained from early-time observations, when the photospheric 
regio ns are partially ionised dPessart & Hillier 2006; Baron et ail 
2007), or from nebular-phase observations during which the pro- 
genitor core, whose properties differ greatly with progenitor main- 
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Table 5. Same as Table |4] but now for the time sequence based on model 
s25el2. 



s25el2 



Age 


^basc.cs 


^phot 


Vphot 


Tphot 


AM phot 


^bol 


rdi 




[10 15 cm] 


[km s" 1 ] 


TK1 


rMni 

L^UJ 


riO 8 Lr^l 


24.0 


2124.66 


1.35 


6524 


8423 


0.040 


12.64 


26.4 


1709.78 


1.46 


6419 


8050 


0.048 


11.87 


29.1 


1356.11 


1.60 


6344 


7815 


0.057 


11.50 


32.0 


1080.90 


1.73 


6274 


7635 


0.067 


11.21 


35.2 


849.14 


1.89 


6203 


7426 


0.080 


10.91 


38.7 


665.03 


2.05 


6121 


7097 


0.098 


10.55 


42.6 


523.98 


2.21 


6009 


6714 


0.128 


10.19 


46.9 


414.59 


2.37 


5858 


6385 


0.179 


9.81 


51.6 


324.49 


2.53 


5668 


6116 


0.263 


9.42 


56.8 


251.96 


2.67 


5437 


5889 


0.405 


9.05 


62.5 


197.03 


2.79 


5166 


5696 


0.639 


8.69 


68.5 


157.38 


2.88 


4866 


5540 


0.997 


8.33 


75.0 


123.25 


2.94 


4535 


5422 


1.518 


7.98 


82.0 


94.47 


2.96 


4183 


5354 


2.243 


7.64 


90.0 


67.84 


2.96 


3803 


5349 


3.244 


7.28 


99.0 


45.06 


2.90 


3385 


5427 


4.480 


6.86 


105.0 


33.80 


2.75 


3028 


5467 


5.473 


6.55 


112.0 


24.62 


2.52 


2605 


5571 


6.457 


5.98 


120.0 


16.11 


2.19 


2117 


5633 


7.317 


5.43 


130.0 


8.67 


1.94 


1723 


5781 


9.138 


4.81 


143.0 


3.30 


1.74 


1410 


5732 


11.628 


3.22 


157.0 


1.16 


1.57 


1154 


5434 


12.923 


1.91 


172.0 


0.84 


1.03 


693 


5270 


13.365 


1.51 


189.0 


0.66 


0.39 


238 


5181 


13.473 


1.30 



sequen ce mass, is revealed dWooslev et all 120021 : | Dessart et al. 
l2010ah . 

The location of the SN II-P photosphere in the outer < 1 Mq 
of the ejecta suggests that a mass buffer of ~ 10 M Q of material sep- 
arates the radiating layer from the birth place of the explosion. As 
long as the photosphere remains in those outer ejecta shells, which 
seems to be during the first half of the plateau phase, radioactive 
decay cannot influence the emergent radiation of SNe II-P. Alto- 
gether, these properties give support to the notion that distance de- 
terminations to SNe II-P using spectroscopy ar e best determined 
from early- rather than late-time observati ons l lBaron et alj|2007l ; 
iDessart et al]|2008l ; iDessart & Hillierll2009h . 

A comparison with simulations of SN 1987 A 
dPessart & H illier 2010) illustrates further the ingredients con- 
trolling the behaviour of the photosphere. At early post-breakout 
times, the rapid recombination of the highly-ionised fast-expanding 
outer ejecta leads to a steep decrease of the photospheric velocity 
(Fig. [4j . In SNe II-pec, this fast decrease is caused in part by the 
relatively small progenitor-star radius, and it slows down only 
when the photosphere reaches the inner ejecta heated by decay 
energy (at ~30d in SN1987A). In SNe II-P, a similar, albeit more 
modest, decrease i n photospheric velocity is seen at early times 
jGezarietal]|2008l) . but it slows down as recombination sets in. 
A rapid decrease kicks in again when the photosphere reaches 
the base of the hydrogen envelope (at JjlOOd in the two present 
simulations). 

Compared to our analyses of SNe II-P 1999em and 2006bp 
( IDessart & H illier 2006]; IDessart et alj|2008l) . characterised by ex- 
plosion ener gies close to but likely larger than the 1.2 B value of 
our models ( Utrobin 20071 ; IDessart et al]l2010al) . the photospheric 









r 




s25el2 i 






sl5el2 * 






=<t- 2 ~~ - 



10 100 1000 

Time since Explosion [d] 



Figure 7. Evolution of the electron-scattering optical depth at the base of 
the ejecta for the simulations based on model sl5el2 (black) and s25el2 
(red). Dots represent the epochs at which our computations are carried out. 
We draw a line of 1 /t 2 slope, which gives the rate of change of the optical 
depth for a homologously expanding gas of non-changing ionisation. The 
faster decline that we obtain in both simulations is caused by ionisation 
changes (i.e., recombination, and primarily of hydrogen), while the step at 
~100d is caused by a composition change as the photosphere enters the 
hydrogen-deficient helium-rich core. 

quantities we obtain with these two sets of simulations are in good 
agreement for the temperature and the velocity. However, our mod- 
els reach larger maximum radii of 2.5-3x10 cm, compared to 
~1.6x 10 cm ( IDessart et alj2008h . We speculate on the origin of 
this discrepancy in §|7] 

As the homologously expanding SN evolves through the 
plateau phase, the ejecta electron-scattering optical depth decreases 
(Fig [7}. It does so initially at constant ionisation, so its rate of de- 
cline follows a 1 /t 2 slope, but as recombination sets in after a few 
weeks, this rate becomes steeper. At about 100 d after explosion, 
the optical depth suddenly declines much faster in both simula- 
tions, as the photosphere enters the hydrogen-deficient helium-rich 
regions of the ejecta, which possess a higher mean-atomic weight 
and hence a smaller mass-absorption coefficient. 

The photosphere vanishes when the total ejecta electron- 
scattering optical depth drops below 2/3 (last entry in Tables [4] 
and |5) and the SN enters the nebular phase. The ejecta are opti- 
cally thin initially in the optical continuum, although not yet in 
the lines. All ejecta mass shells are visible, with emission biased 
towards the denser/hotter inner layers. During the nebular phase, 
the maximum electron density decreases to values on the order of 
10 5 -10 6 cm" 3 , and forbidden lines such as [Ol] 6303-6363A and 
[Call] 7291-7323 A become strong. The nebular ejecta does not 
change its ionisation much with time, and the rate of change of 
the now very low electron-scattering optical depth follows approx- 
imately a l/t 2 slope. For completeness, we show the full evolution 
of the ejecta temperature and electron density over the whole time 
sequence for the simulation based on model sl5el2 in Fig.[6] 



5 SYNTHETIC BOLOMETRIC LUMINOSITY AND 
LIGHT CURVES 

In Fig. [9] we show the bolometric luminosity for each time se- 
quence from 10-20 days until about 1000 days after the explosion. 
As expected for the explosion of RSG progenitors, the initial photo- 
spheric phase is characterised by a quasi-plateau luminosity, with a 
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Figure 8. Left column: Variation of the departure coefficients of Hi (top), Ol (middle), and Fell levels versus Rosseland optical depth for models sl5el2iso at 
10.2 d after shock breakout. Right column: Same as left, but for the model sl5el2iso at 43 d after shock breakout. Non-LTE level populations computed in our 
simulations match closely their LTE counterpart only in regions where the Rosseland mean optical depth is greater than a few. Non-LTE prevails everywhere 
else. 



length on the order of 100-150 d. The longer plateau duration in the 
simulation based on model s25el2 is caused in part by the larger 
production of 56 Ni during the explosion (i.e., 0.163 Mq) compared 
to that in model sl5el2 (0.0866 Mq). The bolometric luminosity 
decreases at the end of the plateau phase and meets the (dashed) 
curve representing the energy-production rate through radioactive 



decay at ~140d (160d) for the simulation based on model sl5el2 
(s25el2)Q 

In our models the fading from the plateau is not abrupt, but 
instead takes about a month, and corresponds to a luminosity de- 
crease by a factor of a few. We find that this prolonged fall-off is 
in part an artefact of the large time step of ~12d (i.e., 10% of 120 
days) in our simulations at such times. Recomputing this phase with 



4 The decay-energy production rate shown in Fig. [9] corresponds to the 
average energy produced per time interval, and thus differs from the instan- 
taneous value at a given post-explosion time. 
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CMFGEN using a finer and constant time increment of 2 d yields a 
faster fading of ~10d but still gives the same fading magnitude 
(Fig.llOt. Our neglect of light-travel-time effects is secondary here, 
since it would only occasion a smearing of the bolometric lumi- 
nosity over a timescale of <ld. Comparison between the evolu- 
tion of the bolometric luminosity and the composition at the photo- 
sphere (Figs.[TJand[9) indicates that the plateau or high-brightness 
phase of the SN corresponds to the epoch during which the photo- 
sphere resides within the hydrogen-rich region of the SN ejecta. As 
soon as the photosphere reaches the outer edge of the helium core, 
the plateau brightness ebbs. The subsequent light-curve evolution 
is qualitatively similar in both models, suggesting little sensitivity 
to the differing helium-core masses (4.27 Mq in model sl5el2 and 
8. 19 M Q in model s25el2). 

Bolometric light curves for four SN II-P (S N 1999 e m, SN 
2003gd, SN 2004et, SN 2005cs) are presented by ISmarttl J2009T) . 
All four show a drop of around 1 magnitude close to 120d from 
maximum, and for one of them (SN 2005cs) the bolometric lumi- 
nosity changes by almost a factor of 2 on a time scale of ~ 10 days, 
while for SN 1999em the change is more gradual. As the end of the 
plateau phase is determined by the transition of the photosphere 
from the hydrogen rich core to the helium core, the steepness of the 
change will depend on any mixing that has occurred between these 
layers, as well any departure from spherical symmetry. 

The bolometric luminosity for the simulation based on model 
s25el2 is greater at all times than for the simulation based on model 
sl5el2. The larger plateau brightness stems from the larger progen- 
itor radius in model s25el2, while the larger nebular flux is caused 
by the larger amount of 56 Ni produced by the explosion in model 
s25el2. The magnitude of the fading at the end of the plateau is 
thus a function of both the progenitor size and the amount of 56 Ni 
produced. As the progenitor radius and mass of nucleosynthesised 
56 Ni are unlikely to be fundamentally related (because of the in- 
fluence of abundance variations, rotation, binary evolution, and the 
complexities of the explosion mechanism), there should be no clear 
correlation between plateau and nebular brightness, and thus the 
observed post-plateau fading should show a large scatter. The ap- 
proximate treatment of radioactive decay in this study cannot alter 
this result since no 7-rays are predicted to escape up to >500 d. At 
later times, our neglect of 7-ray photon escape will lead to an over- 
estimate of the bolometric luminosity shown here, by up to a factor 
of two at 1000 d. 

The bolometric luminosity is not directly observable because 
not all spectral bands are accessible. Moreover, filter bandpasses 
are separated by gaps, sometimes sizeable, that miss regions of 
non-zero flux (and potentially strong lines). Regions of intense at- 
mospheric absorption prevent the observation of certain spectral 
regions altogether, as in some parts of the near-infrared (near-IR). 
Hence, in Fig. QTJ we show multi-band synthetic photometry for 
our two simulations based on models sl5el2 and s25el2. These 
are obtained by integrating our non-LTE time-dependent synthetic 
spectra at eac h epoch with the filter transmission functions of 
lLandoltl d 1992b for the optical bands, and with those of Co hen et"aTJ 
1 20031) for the near-IR bands. 

While the bolometric luminosity shows a steady and slow de- 
cline at most epochs, the multi-band light-curves show a contrast- 
ing behaviour affected by ejecta cooling. In particular, the photo- 
spheric cooling after shock breakout leads to a continuous decrease 
of the flux at short wavelengths, as the peak of the distribution 
sweep s through the far-U V, to the UV, and approaches the optical 
bands dGezari et alj|2008l) . This translates here into a "bell" shape 
for the B— band light curve, but a plateau followed by a steady fad- 
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Figure 9. Main panel: Log-log plot of the evolution of the bolometric lumi- 
nosity for our simulations based on model sl5el2 (black; model contains 
0.0866 M of 56 Ni initially) and s25el2 (red; model contains 0.163 M of 
56 Ni initially). Dots indicate the time of each computation. We also over- 
plot as a dashed line the mean 7-ray energy released per time interval, com- 
puted independently for each model using our Monte Carlo code (see Hillier 
& Dessart in prep.), and assumed in this set of calculations to be deposited 
locally as a heat source (we do not treat non-thermal ionisation/excitation). 
Note that due to computation over different time intervals, the solid/dashed 
curves diverge somewhat at very late times, when time increments are large 
and on the order of 100 days. Insert: Same as main panel, but now zooming- 
in on the plateau phase and the onset of the nebular phase. We now use a 
linear scale for both axes. 
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Figure 10. Same as Fig. [5] for models sl5el2 and s25el2, but now com- 
paring results from CMFGEN when we adopt a time increment of 10% of 
the current time (black; this is our standard choice already shown in Fig. [9] 
which yields large time steps at late times), or a constant time increment of 
2 to 5 days (blue curves). For each curve, the dots give the times at which 
the CMFGEN calculations were carried out. Using a fixed and finer time in- 
crement yields a different result at the end of the plateau phase, which is 
better resolved and found to last about 10 d in both models. 



ing in the [/—band. In contrast, the VRIJHK magnitudes all rise 
at all times prior to about 50 days after explosion, before turning 
over and declining. Both simulations based on models sl5el2 and 
s25el2 show a similar behaviour, the latter being visually brighter 
at all photospheric and nebular epochs. As discussed previously, 
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Figure 11. Top: Synthetic optical and near-IR light curves for the simu- 
lation based on model sl5el2 and obtained by integrating our non-LTE 
ti me-dependent sp ectr a at each epoch over the filter-transmission functions 
of lLandoltl 1 19921) and lCohen et alj (2003h . Bottom: Same as top, but now 
for the simulation based on model s25el2. Dots indicate the time of each 
computation. 



this contrast stems from the larger progenitor radius, and the larger 
amount of 56 Ni produced in model s25el2. 

Besides differences in progenitor properties, we also find 
that magnitudes/colors are affected by the completeness of the 
model atom employed for the calculations. Model sl5el2iso is 
identical in ejecta structure to that of model sl5el2 (same ra- 
dius/velocity/density/temperature distribution), but differs in com- 
position due to the treatment of under-abundant species, ignored 
in model sl5el2. We find stronger line blanketing in simulations 
based on model sl5el2iso at the recombination epoch, in partic- 
ular associated with the treatment of Sc II, Cr II, Fe I, and Ni II, 
which leads to a reduction of the already faint (7/3-band magni- 
tudes b y <CQ.5mag. As already realised for the modelling of SN 
1987A jPessart & Hillienl2010h . including all species, and in par- 
ticular all under-abundant intermediate-mass-elements/metals up to 
Fe/Co/Ni, and as many ionisation stages as possible is critical for 
obtaining converged results. In SNe II modelling, however, it is not 
critical to include a huge number of levels, as we find that increas- 
ing the number of levels by a factor of a few from the choice given 
in Table[3]does not alter our results significantly. 



6 SYNTHETIC SPECTRA 
6.1 Global evolution 

We now turn to a presentation of the synthetic spectra for our two 
simulations. Because the general spectroscopic evolution from the 
photospheric phase to the nebular phase is similar for both se- 
quences, we show in Fig. [12] a montage for one simulation only, 
that based on model sl5el2. Our continuous modelling is permit- 
ted by the treatment of the entire ejecta at all times. The inner-grid 
boundary is at the inner edge of the innermost ejecta shell; any ma- 
terial between that shell and the neutron star fell back and is thus 
not part of the ejecta by definition. In the full solution of the ra- 
diative transfer with CMFGEN, we adopt two inner-boundary types. 
The inner boundary is treated as diffusive (with an imposed inner- 
boundary luminosity of Lq/10 6 ) as long as the base optical depth 
is greater than ~10. Afterwards, we switch to a zero-flux nebular- 
like inner-boundary. For the computation of spectra, we assume a 
hollow core. With such a full-ejecta approach, we can continuously 
model the ejecta as they evolve from optically thick to optically 
thin. 

Figure[T2]shows the evolution of a spectral-energy distribution 
that initially resembles that of a thermal emitter, with a decreasing 
continuum flux at 100-150 d after explosion when the ejecta be- 
come optically thin. As time progresses in the nebular phase, the P- 
Cygni profile morphology of lines transitions to that of pure emis- 
sion, with a box-like structure typical of optically-thin line forma- 
tion f rom a hollow region expanding at near-constant speed ( ICastod 
1970). Simultaneously, the continuum flux vanishes. 

During the SN evolution, line features evolve considerably, re- 
flecting the change in composition, ionisation/temperature and den- 
sity. At early times, the photospheric layers are hot and ionised and 
H I Balmer lines and He 1 5878A are present. As recombination sets 
in, line blanketing caused by once-ionised intermediate-mass ele- 
ments and metal species appears at all wavelengths smaller than 
~5000 A (Fig.l 13t. This entire photospheric-phase evolution takes 
place at essentially constant composition, so that these spectral 
changes are governed by ionisation changes. As the SN progresses 
into the nebular phase, lines forming through recombination (fol- 
lowing photoionisation) weaken because of the fading radiation 
field, and in particular the quasi absence of ionising photons, either 
in the Lyman or in the Balmer continuum (in reality, this is exagger- 
ated by our neglect of non-thermal ionisation by high-energy elec- 
trons). In contrast, forbidden lines, which form at low density by 
radiative de-excitation of collisionally-excited levels, strengthen on 
top of an ever-fading continuum. The two most important forbid- 
den lines that appear in essentially all nebular-phase SNe II spectra 
are the doublets [O I] 6303-6363 A and [Ca II] 7291-7323 A. 



6.2 Nebular phase 

Due to their faintness, and unimportance for cosmological studies, 
there have been relatively few nebular studies of SN. The most ex- 
tensive nebular studies are for SN1987A. For example. [Xu et all 
dl99lh studied the Compto nization of gamma rays by cold elec- 
trons. lxu & McCravl ( fl99ll) studied how non-thermal electrons are 
produced and d owngraded. IXu et al. (l992) studied the formation 
of the Hi lines, lli & McCravl 1 1995 ) stu died the form ation of the 
He I emission lines, while iLi & McCravl d 1 995L 1 1 993l> studied the 
formation of the Call emission lines. Additi onal extensive stud- 
ies of t he nebular phase in SN1987A are by iKozma & Franssonl 
dl998albb who discussed in detail the temperature and ionisation 
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Figure 12. Montage of synthetic optical spectra covering epochs between 10 and 1000 d after explosion for our simulation based on model sl5el2. The 
selected times, indicated by the label on the right, are such that the time interval between two consecutive spectra is constant in the log. For better visibility, 
each spectrum is first normalised to the peak flux in the region 6000-8000A, which introduces a variable scaling between epochs, and then shifted vertically 
for better visibility. During the photospheric phase, lines features are associated with Hi Balmer lines, resonance lines like the triplet of Call at ~8500 A and 
forests of metals lines (e.g., Fell and Till) causing severe line blanketing at wavelengths smaller than ~5000 A. During the nebular phase, forbidden lines like 
[Ol] 6303-6363 A and [Call] 7291-7323 A progressively strengthen relative to recombination lines and the fading continuum. 
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Figure 13. Illustration showing bound-bound contributions of selected species (bottom; fluxes are normalised to the continuum and shifted vertically) to the 
full synthetic spectrum (red curve, top panel) obtained in the simulation based on model sl5el2iso at a post-explosion time of 75.8 d. The full and continuum- 
only (blue curve) synthetic spectra have been reddened (E(B — V) = 0.1) and distance scaled (D = 11.5 Mpc; an additionnal scaling by a factor of 2.3 is 
added to compens ate for the higher br ightness of the model), and in the top panel are compared for illustration to the Keck observations of SN1999em on 5 
Dec. 1999 (black: iLeonard et al.l20o1 
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Figure 14. Illustration of optical (left column) and near-IR (right column) synthetic spectra obtained for the simulation based on model sl5el2 at post- 
explosion times of 302 d (top row) and 1040 d (bottom row). The full synthetic spectrum is shown as a black curve, while overplotted coloured lines refer to 
the flux obtained when including bound-bound transitions exclusively of the species labelled on the right of each panel. We only display the results for those 
species that matter to the emergent flux, hence the selection varies between panels (i.e., for different times and spectral ranges). At late times, only forbidden 
lines, forming in the inner SN-ejecta regions corresponding to the former progenitor helium core, contribute to the SN light. The continuum flux is a few orders 
of magnitude smaller than the full synthetic flux due to lines. 



structure, and the formation of the most important emission lines. 
A major result of the above studies is the need for mixing of hydro- 
gen down to 70 km s -1 , with oxygen and h elium present at even 
lower velocities. iKozma & FranssorJ ql998bh analyzed the O I lines 
in SN 1987A and found that the bulk of the oxygen was located be- 
tween 400 and 2 000 km s -1 , a result in broad agreement with the 
earlier studies of iLi & McCravl dl992l) . 

During the nebular phase the continuum is optically thin and 
line formation is simplified. Line profiles during the nebular phase 
provide information about the distribution of material as a func- 
tion of velocity, and its geometrical distribution. For example, an 
isolated emission line arising from a spherical hollow shell will ex- 
hibit a flat top, with its minimum width determine d by the velocity 
at the inner edge of the shell (e.g.. lMcCravlll993h . As the Sobolev 
escape probability is independent of angle in a Hubble flow (see 
expressions of ICastorl 1970), this result holds independent of the 
line's optical depth (neglecting continuum absorption). 

During the nebular phase, one probes the entire 
optically-thin ejecta, although emission is biased towards the 
denser/slower/hotter inner regions where unstable isotopes are 
naturally more abundant. This inner region contains few tens of 
percent of the total ejecta mass, and <10% of the total ejecta 
kinetic energy in our models (Fig. QJ. The faintness of all SNe at 
such epochs makes them difficult to observe at late times. However, 



because the nebular phase reveals the core of SNe II-P progenitors, 
it is a key source of information on the explosion, as well as the 
core properties of the progenitor star. Despite this importance, 
nebular-phase spectra of SNe II-P are scarce and th us not much 
has be en extracted yet from this important phase foessart et al.l 
l2010al) . Dedicated campaigns to secure nebular spectra of SNe II-P 
are needed to improve upon this matter. 

We show in Fig. [14] the optical and near-IR synthetic spectra 
obtained for the simulation based on the model sl5el2 at 302 and 
1040 d after explosion. While the black curve gives the full spec- 
trum, coloured curves give the bound-bound contributions from in- 
dividual species. At such late times, not only 56 Ni, but also 56 Co, 
have decayed and the metal mixture is dominated by 56 Fe. At 
302 d, the synthetic spectrum is dominated by a forest of Ti II and 
Fe II lines in the region blue-ward of 5500 A. Further to the red, 
fewer Fell lines are present (though in the simulation based on 
model sl5el2iso, lines of Fel are predicted), and three main fea- 
tures dominate the spectral landscape, the forbidden doublet lines 
of [O I] 6303-6363 A and [Can] 7291-7323 A, and the Can triplet 
at 8498, 8542, 8662 A. As time proceeds, the decreasing ejecta 
electron/ion density favours the forbidden and "resonance" tran- 
sitions and so the features associated with [Ol], [Call], and Call 
strengthen while the recombination transitions associated with Fell 
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weaken in the region 4000-5500 A. At very late times, only [O I] 
and [Fe II] persist. 

In our simulations we do not predict Hi Balmer-line emis- 
sion at nebular times. The Balmer lines fade abruptly when the 
ejecta become optically thin, immediately at the end of the plateau 
phase, a result of the sudden disappearance of Balmer-continuum 
photons and the de-excitation of H I. The treatment of non-thermal 
ionisation/excitation would counteract the fading of the Hi Balmer 
recombination lines at the end of the plateau phase, since it 
is precisely at that time that high-energy electrons can impact 
level populations and the ionisa tion state of the gas in SNe II 
dKozma & Franssonl 19921.1 1998allbh . We suspect this neglect is less 
critical for forbidden lines which form through transitions between 
levels that have a low excitation potential, and hence can be excited 
by thermal electrons. 

The rapid disappearance of the Hi lines at the end of the 
plateau phase indicates that non-thermal processes must be crucial 
for the hydrogen ionisation, even at ~ 150 days. This is in a pparen t 
contradiction with the conclusions of Koz ma & Franssonl (1998b) 
who indicate that in SN1987A, H ionisation is determined by pho- 
toionisation from the n=2 level up to 500 d. A likely explanation 
is that non-thermal excitation from the ground state, together with 
the large optical of Lya, helps to maintain (set) the n=2 popula- 
tion from which ionisation occurs. The viability of this option can 
be verified by solving th e two-level atom plus con tinuum problem 
(see, e.g., equation 39 of iKozma & Fransson 1993). 

In model sl5el2, and at a representative post-explosion time 
of 332 d, we find that the doublet lines [O I] 6303-6363 A and 
[Ca II] 7291-7323 A form primarily in the inner denser ejecta. This 
is illustrated by the similar and relatively modest broadening, of 
Doppler origin, of the corresponding individual line components. 
And indeed, the full-widfh-at-half-maximum of ~2500km s _1 
corresponds quite closely to the velocity of ~1500km s _1 of the 
outer helium-core (see Fig.[T] but also Fig. [2j. 

In contrast, the strong permitted transitions of Call at 8498, 
8542, 8662 A show a hybrid formation, with contributions from 
both the inner and outer ejecta. This gives a very peculiar shape to 
the Ca II 8662 A line, with a relatively narrow emission part associ- 
ated with the ejecta region corresponding to the progenitor helium- 
core, and a broad emission part associated with the hydrogen-rich 
region of the SN ejecta (Fig. I15t . The second, broad, component 
suggests that emission comes from the entire volume of the ejecta, 
all the way to the outermost mass shells (see Fig. QJ. Excluding 
Ca II 8662 A line (discussed below), simulated emission line pro- 
files show that the bulk of their emission at nebular times comes 
from the same region of space, which corresponds to the progeni- 
tor helium core. 

The Grotrian diagram for Call is rather simple and it is not 
initially obvious why Call 8662 A should show such a different 
line profile from the other lines. A hint comes from comparing the 
Call triplet spectrum to that obtained by summing the individual 
profiles. As shown in Fig. [T6] there is a lot more flux in the two 
bluest lines than in the full spectrum. This suggests that the photons 
emitted in Call at 8498 and 8542 are scattered in the Call 8662 A 
elsewhere in the envelope. This is reasonable since the triplet transi- 
tions are optically thick, and the velocity shift between neighbour- 
ing lines is only ~ 1500 km s _1 and ~ 4200 km s - . Examination 
of the Call triplet profile on day 198 for SN1987A also shows a 
complex profile, wi th the 8662A line being much broader th an the 
8498/8542 A blend (Phillips et alJl99d : lLi & McCravll 1993b . 

Our theoretical spectra also reveal a wealth of weak emission 
lines which blend with many of the stronger features. This was pre- 
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Figure 15. Illustration of the synthetic line flux for the simulation based 
on model sl5el2 at 332 d after explosion, and associated with individual 
transitions giving rise to the forbidden lines of [O I] and [Call], as well as 
the "resonance" transitions of Ca II appearing around 8500 A. All fluxes are 
normalised to the peak flux of the Ol 6303 line. A colour coding is used 
to differentiate individual line transitions. Notice the similar and relatively 
narrow width of most of these lines, reflecting the expansion speed of the 
helium-core in which they form, while resonance transitions of Call form 
not only in that same inner region, but also in the outer hydrogen-rich ejecta. 
Individual profiles were computed ignoring interactions with overlapping 
lines. 
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Figure 16. A comparison of the theoretical spectrum in the neighbourhood 
of the Call triplet ("FULL"; green), and that obtained by summing line pro- 
files for individual members of the triplet ("SUM"; red). Individual com- 
ponents are drawn as thin coloured lines. The "summed" spectrum is very 
different from the true spectrum, particularly for the 8498 and 8542 A lines 
which contain too much flux in the "summed" spectrum. As explained in 
the text, photons in 8498 and 8542 A are absorbed and re-emitted in the 
8662 A line. 



viously highlighted by IKozma & Franssonl dl998bl) for O I. They 
noted that the ratio of O l A6300 to A6363 was significantly less 
than the theoretical value of 3., which they attribute to a blend of 
Fe I 6361 A. In our calculations, we find that this [Ol] region over- 
laps with numerous Fel and Fell lines (Fig.|17t, and t hus may com- 
plicate the determination of the oxygen-mass ejected (Dessart et al.l 
l2010al) . 
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Figure 17. Full synthetic spectrum (black) for model sl5el2iso at 306 d 
after explosion and showing the [Ol] 6303-6363A region. To illustrate the 
presence of line overlap, we overplot the synthetic spectrum obtained by ac- 
counting for bound-bound transitions of Ol (blue), Fel (green), or Fell (red) 
exclusively. Line overlap and associated line-optical depth effects compro- 
mise any eye-ball estimate of the flux contribution associated with Ol. 

6.3 Comparison between simulations based on the sl5 and 
s25 models 

We have so far focused on the spectroscopic properties obtained 
for one time sequence. We now compare the spectroscopic differ- 
ences between the simulations based on models sl5el2 and s25el2, 
which both have the same explosion energy, were evolved under the 
same conditions (no rotation, solar-metallicity environment), but 
started from stars with main-sequence masses of 15 and 25 Mq. In 
§[2] we discussed the key differences between the two pre-SN stars, 
which are the larger progenitor radius in model s25el2, the more 
pronounced He/N enhancement in the envelope of model s25el2, 
and the larger helium-core mass in model s25el2. Indeed, model 
s25el2 ejects 6.29 Mq of helium-core embedded material, while 
sl5el2 ejects only 2.44 Mq. Here, we discuss the radiative signa- 
tures resulting from such differences in progenitor and explosion 
properties. We show comparative spectra for the two simulations at 
the same post-explosion time during the photospheric and nebular 
phases in Fig.l 1 81 

As discussed above, the larger progenitor radius for a given 
explosion energy gives bluer colours and a brighter bolometric dis- 
play (larger flux) during the plateau phase, while the large 56 Ni 
yields in model s25el2 give a higher nebular luminosity (at a given 
post-explosion time). During the recombination epoch, when the 
simulations show similar colors, contemporaneous spectra obtained 
for the two simulations show a remarkable similarity. The differ- 
ence in composition, imprinted in H/He/C/N/O is largely masked 
and photospheric-phase spectroscopy is unable to lift this apparent 
degeneracy. The few C I, N I and O I lines appearing in the red part 
of the optical are too weak to allow a reliable and meaningful as- 
sessment of CNO abundances. Past about 15 d after explosion, Hel 
lines disappear, and the weak sensitivity of the H I Balmer lines to 
the H/He abundance ratio compromises even this determination. 

Quite paradoxically, most inferences on SNe II-P are based 
on observations during the photospheric phase, which, as we have 
illustrated here, boasts very degenerate properties that prevent a 
clear assessment of the progenitor properties. Having recourse to 
nebular-phase spectra can help lift this degeneracy in a number of 
ways. Because the helium-core mass increases with main-sequence 
mass, the relative fraction of the ejecta that it represents grows with 



main-sequence mass. iDessart et alj J2010ah demonstrated that for 
a given explosion energy, the ejection velocity of core-embedded 
oxygen-rich material is an increasing function of main-sequence 
mass. The width of [O I] 6303-6363 A can be used as a diagnostic 
of this differential expansion speed. Its strength can then be used to 
constrain the total amount of core-embedded oxygen ejected. To- 
gether with the constraint on the helium-core mass (which is set by 
the main-sequence mass), one may have a means to constrain the 
mass of the compact remnant formed. 

We plan a forthcoming parameter study using non-LTE time- 
dependent simulations based on models of SNe II-P with different 
main-sequence mass, explosion energies, and 56 Ni. This will re- 
veal a wealth of information, and document what SNe II-P ejecta 
properties condition the [O I] strength and width. Based on the two 
simulations presented here, we find that the doublet line [O I] 6303- 
6363 A is both broader and stronger in the simulation based on 
the higher-mass progenitor, and this results from differences in 
the pre-SN stars rather than in the explosion properties (Fig. [T9j 
note that the enforcement of homology alters somewhat the inner- 
ejecta kinematics). The ~ 30% larger core-embedded oxygen ejec- 
tion speed in model s25el2 is reflected by the greater width of 
each component of the [Ol] doublet line, while the larger core- 
embedded oxygen ejected mass in model s25el2 (3.32 compared to 
O.8I6M0) is reflected by the higher line flux at all nebular epochs 
(this is also influenced by the 56 Ni mass). Clumping and subtleties 
associated with radioactive heating will complicate a quantitative 
assessment here, but this differential reasoning at least opens a 
route towards constraining in a systematic fashion the properties 
of the inner ejecta of SNe II-P. 

6.4 Comparison between simulations based on sl5el2 and 
sl5el2iso models 

As detailed previously in §f2] models sl5el2 and sl5el2iso differ 
in their composition, and our simulations based on these simula- 
tions reflect that distinction, as well as the additional treatment of 
ions in the model atom (see Table [3} . In our simulation based on 
model sl5el2, we also neglected some species for which no asso- 
ciated line feature was ever seen, such as neon or argon. In that 
case, we applied a scaling on the rest of the elements to normalise 
the sum of mass fractions to unity - this scaling modifies by ~10% 
the abundances in the helium core. The difference in composi- 
tion between the two models contains two components. First, some 
species are included in model sl5el2iso, like Sc or Ba, which are 
missing in model sl5el2. Second, model sl5el2iso gives the distri- 
bution of species versus depth. It shows, for example, an increase in 
sodium abundance towards the core, while in model sl5el2 sodium 
is taken at its solar-metallicity value. 

So, when comparing synthetic spectra from our simulations 
based on models sl5el2 and sl5el2iso, both composition and 
model atom differences intervene. We show in Fig. [20] the differ- 
ences at the same post-explosion time of 75 d for each simulation. 
There is ~5% absolute flux level difference between the two, with 
significant differences appearing in the UV, and essentially none 
beyond <;1 /im. The enhanced UV flux in the simulations based on 
model sl5el2iso may stem from the reduced metal abundances, in 
favour of the treatment of species like Ar or Ne which have lit- 
tle blocking power, while the reduced U B-band flux is caused by 
the treatment of additional blanketing sources (e.g., Fel or Nill; 
see Fig. [T3] for an illustration of individual contributions). Lines 
of Sell, whose abundance is zero in model sl5el2, are predicted 
in the simulation based on model s!5el2iso. Similarly, the NalD 
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Figure 18. Comparison between optical synthetic spectra for the simulations based on model sl5el2 (red) and s25el2 (black) and shown at post-explosion 
times of 25 d (top left), HOd (top right), 300 d (bottom left), and 1000 d (bottom right). Only a distance scaling is performed on these synthetic spectra, 
adopting a value of lOMpc. 
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Figure 19. Same as Fig[l8] but now showing the [O I] i 

doublet line is stronger in that simulation, which reflects the in- 
creasing sodium mass fraction in model sl5el2iso with depth. Dur- 
ing the nebular phase, differences are significant between the two 
simulations (not shown), but they are primarily associated with the 
treatment of Fe I. 

Although the spectral evolution computed for simulations 
based on models sl5el2 and sl5el2iso are essentially the same, the 
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1-6363A doublet region at 300 d (left) and 1000 d (right). 

model with a more complete set of metals and model atom yields 
stronger line blanketing (which translates into a weaker flux in the 
UB bands) as well as the prediction of specific lines associated 
with, for example, Sell (Fig. 113) , 
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Figure 21. Left: Comparison between the observed m ulti-color light-cur ve (diamonds) of SN1999em iLeonard et ail 2002) and the con'esponding reddened 
(E(B - V) =0.1) a nd distance sca l ed (D= l 1.5 Mpc; IBessart & Hillieill2006h magnitudes from the sl5el2 model (dots). A further vertical shift is applied 
to avoid overlap, as in[Leonard et al. 12002). Note the excellent agreement at day 9 for each band (see spectral comparison at that epoch in Fig. [22), but the 
subsequent divergence between observed and synthetic magnitudes. This suggests a significant mismatch between the structure of the RSG progenitor model 
and/or explosion energy with that corresponding to the progenitor of 99em. This may stem from the mismatch in density distribution, which is in evidence in 
Fig. 1221 Right: Same as left, but now for the time sequence based on model s25el2 (we now scale the synthetic magnitudes adopting a SN distance of 17 Mpc). 



7 COMPARISON WITH OBSERVATIONS 



7.1 Bolometric luminosity and light curves 



Observed l/-band ligh t curves of a lar ge sample of SNe II-P 
have been presented in lHamuvl d2003h or Poznanski et a D d2009l) . 
The se suggest that desp ite the expected variatio ns in progeni- 
tor dLevesque et al]|2005t) and explosion properties dBurrows et alj 
2007b. a; Marek & Janka 2009), their properties are quite uniform 
(those wishing to standardise such SNe for distance determinations 
will find instead that there is a sizeable scatter). Quite generally, 
about 10-20 days after the inferred time of explosion, the visual 
brightness levels off or fades slowly (corresponding to a typical 
luminosity of ~1.2x 10 42 erg). When reaching 75-120 d after ex- 
plosion, the SN fades by 0.8-1.2 dex in luminosity as it transitions 
on a ti mescale of a few weeks fr om the plateau phase to the nebular 
phase dBersten & Hamuy||2009l) . 

Our simulations reveal plateau-phase bolometric luminosities 



(Fig. [9} that are a factor of two too large compared to obser- 
vations. Indeed, the absolute IZ-band light curves we obtain are 
brighter than those we derived from detailed modelling of SN e 
1999em and 2006bp dDessart & Hillierll2006l: fPessart et alj|2008h . 
To illustrate this discrepancy, we show a comparison of our red- 
dened and d istance-scaled synth etic light curves to observations of 
SN1999em dLeonard et al.l20 02) in Fig. [21] Interestingly, the same 
discrepan cy applies to the simulatio ns of SNe II-P light curves pre- 
sented in iKasen & Wooslevl d2009l) . In contrast, the simulation of 
lEastman et al.l dl994l) . which is based on a 15 Mq progenitor star 
characterised by a smaller radius of 480 Kq, predicts a represen- 
tative plateau-phase bolometric-luminosity of ~4x10 8 Lq and a 
peak V-band magnitude of about —17. Hence, we surmise that the 
main problem with the progenitors used in recent SNe II-P simu- 
lations (including this work) is the progenitor radius, which is sys- 
tematically too large, and perhaps the helium mass fraction in the 
hydrogen envelope, which is too low. 

Our simulations have an explosion energy of 1.2Banda 56 Ni 
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Figure 22. Left: Comparison between the observations of SN1999em on the 3rd of November 1999 <Hamuvetall200lh with the reddened and distance-scaled 
synthetic spectrum obta ined from the sl5el2iso model on day 9.3 evolved with the time-dependent solver. Right: Same as left, but now using the model of 
Dessa rt & Hillie] 1 2006) for the same date, which is computed with a steady-state solver (which also accounts for v/c terms appearing in the radiative-transfer 
equation). The difference between the two models is caused by the steeper density distribution in the hydrodynamical input sl5el2iso, corresponding to 
N p ~ 25 in the photospheric region at the time, compared to N p =10 used in the steady-state calculation. This affects significantly the flux in recombination 
lines, as well as the evolution of the light curve and colours, and suggests that the hydro input is not fully adequate. 



mass that may not fit specific SNe II-P but is representative of the 
standard and thus should not be the cause of this overall discrep- 
ancy. Final ly, in the simulations o f SN1987A performed with the 
same code (Dessart & Hilliej2010h . we did reproduce the observed 
SN1987A light curve to within 10% in flux at all times between 0.3 
and 21 d. This agreement suggests there is a physical discrepancy 
with the pre-SN RSG stars we employ at present. We speculate that 
RSG progenitors with a surface radius intermediate between the 
blue-supergiant progenitor of SN1987A (~50Rq) and those of the 
RSG progenitors of models sl5el2/s25el2 (810/1350 R ) would 
correspondingly reveal light curves that are intermediate between 
those produced from our simulations of SN1987A and those from 
models sl5el2/s25el2. 

A faster increase in (or a globally-larger) helium mass fraction 
in the hydrogen-rich envelope would likely remedy the sustained 
visual brightening obtained with our models during the plateau 
phase. Indeed, in this context, the progression of the recombination 
wave is faster, leading to a faster decline of the photospheric ra- 
dius than presently predicted (Fig. [3j. In ou r non-LTE steady -state 
simulations for SNe 2005cs and 2006bp jPessart et aT]|2008h . we 
employed a chemical mixture with dep leted H/C/O and e nhanced 
He/N (a similar trend was obtained bv lBaron et alj|2007h . resem- 
bling closely the surface composition of model s25el2, but not that 
of model sl5el2. As discussed earlier, there is evidence of mix- 
ing in RSG stars wh ich points in this direction of enhanced He/N 
dLancon et al.l 120071) . perhaps as a result of stellar rotation. This 
determination of H/He/C/N/O abundances is best done with early- 
time observations of SNe II-P, when the photospheric layers are 
fully or partially ionised. This issue is important to resolve, as it 
has ramifications both for the modelling of SN light and for the 
understanding of mixing processes in massive stars prior to the ex- 
plosion. 

As in our simulations of SN1987A, there are residual prob- 
lems with the UV range and [/B-bands. The flux that falls within 
these regions is typically smaller, or much smaller, that the vi- 
sual flux. This region is strongly affected by line blanketing, which 
effectively absorbs the relatively small flux appearing blue-ward 
of <^4000A at the thermalisation depth. This line blanketing is 
caused by under-abundant metal species, and may also be affected 



by the floor temperature we employ in our simulations. For exam- 
ple, the additional contributions from Sell, Cr II, Fe I and Ni II in 
simulations based on model sl5el2iso led to a steeper decline and 
<0.5 mag-fainter U B-band magnitudes than obtained with model 
sl5el2. This sensitivity suggests that the discrepancy is really con- 
nected to completeness of the model atoms, although tests have 
not revealed any obvious candidate species, or a problem with the 
adopted model atoms. As we neglect many neutral species we will 
underestimate the line blanketing, particularly at late time. Exam- 
ination of the 17-band light curve presented by Kasen & Wooslevl 
(2009) for their standard II-P model shows less of a discrepancy 
with the U — V color, further hinting to a possible problem with 
our adopted atomic data. 

The metallicity at the site of the explosion could also play a 
role (it may not be solar, as assumed here), as could the H/He abun- 
dance ratio. The H/He abundance ratio modifies the way radiation 
thermalises so that for enhanced helium abundance, the UV flux 
drops. A higher helium mass fraction in the hydrogen-envelope of 
the RSG progenitors would thus lead to a faster visual-brightness 
decline and a fainter UV flux, both improving the match to standard 
observations of SNe II-P. 



7.2 Photospheric- and nebular-phase spectra 

As described above, the absolute brightness and spectral evolution 
predicted with the simulations performed in this paper differ quan- 
titatively from that of standard SNe II-P. However, modulo an abso- 
lute scaling to get the fluxes in correspondence, and allowing for an 
adjustment in time, we can compare our synthetic spectra to indi- 
vidual observations. Here, we use the spectroscopic observations of 
SN1999em, for which we als o have tailored simulati ons computed 
with a steady-state approach jPessart & Hillierl2006h . 

Our simulation at 9.3 d and based on mode l s!5el2iso matches 
the SN 1999em observations on 3 Nov. 1999 dHamuv et ai]|200lh 
both in relative and in absolute flux in the optical (Fig. 1221 . Sim- 
ilarly, the e quivalent steady-state m odel computed for that date 
(right panel; iDessart & Hillierll2006h shows the same color and re- 
produces the overall flux very well. However, the equivalent steady- 
state model reproduces line profile strengths much better. Time de- 
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pendence might play a role here, but more likely, the mass-density 
profile in the outer ejecta of model sl5el2iso is the culprit. 

At very early post-explosion times, a density exponent of 
15-50 is req uired to reproduce the observations of SNe 2005cs 
and 2006bp jPessart et alj|2008l) , but at 10 days after explosion, 
this density exponent was found to be on the order of 10 in 
the photospheric region for these SNe, as well as for SN1999em 
jPessart & H illier 2006). In our simulations (sl5el2/sl5el2iso), 
the photosphere resides from day 10 to day 25 in the outer ejecta 
layers where the density exponent is approximately 25 (!). This 
then seems the most probable reason for the underestimated line 
strength (see discussion of the eff ect of the density expone nt on 
Balmer-line profile morphology in iDessart & Hiliieill2005bh . The 
progenitor star may have a different structure. For example, the 
density exponent for the corresponding regions (velocity in the 
range 7000-9000 km s" 1 ) in our 87A mo del ejecta, produced from 
the explosion of a 50 R BSG star, is 8 jPessart & Hillierll2010h . 
Furthermore, simulations of SNe II-P ejecta tend to have low res- 
olution in the outermost mass shells, which questions the accuracy 
of these layers and the corresponding early-time SN photospheric 
properties. 

As the ejecta recombine and the photosphere recedes into 
deeper layers, the line strengths increase. At such times, we re- 
produce the observations of SN1999em extremely well (upper two 
panels of Fig. 123 b- For the observations of SN1999em on 5 Dec 
1999, we also illustrate the various bound-bound transitions that 
contribute to the observed line features in Fig.[T3] We also include 
here a comparison of th e nebular-phase sp ectrum of SN1999em 
taken on 26 Sept. 2000 dLeonard et al. I l2002h with the synthetic 
spectrum at 306 d from our simulation based on model sl5el2iso 
(synthetic spectra are scaled to correct for the magnitude mismatch 
illustrated in Fig. l2"Tt . Reca ll that SN199 9em is thought to have 
produced 0.036 M Q of 56 Ni JUtrobinl 20071) while model sl5el2iso 
contains initially 0.086 Mq. The agreement is good, especially 
since no parameter is adjusted in model sl5el2iso to match the 
observations. The only, and sizeable discrepancy, is the complete 
absence of Ha in our synthetic spectrum, while a strong flux is 
observed instead. We also note that we predict a sudden disappear- 
ance of the continuum flux at the end of the plateau phase, while 
observations suggest a more gradual fading of the continuum flux 
relative to lines. Such discrepancies stem entirely from our neglect 
of non-thermal excitation/ionisation. Remarkably, it may be that 
this additional process does not affect other species mucrQ In any 
case, this physical ingredient is essential for quantitative modelling 
of the nebular phase of SNe II and will be incorporated in CMFGEN 
in the near future. 



8 DISCUSSION AND CONCLUSIONS 

We have presented the first non-LTE time-dependent radiative- 
transfer simulations of SNe II-P covering from the early pho- 
tospheric phase until >1000d after explosion. Our simulations 
start from hydrodynamical models of the explosion of two RSG 
stars, evolved from non-rotating solar-metallicity stars with main- 
sequence masses of 15 and 25 M© . The main assets of our approach 
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Figure 23. Top: Comparison between the reddened synthetic spectrum 
from our simulations sl5el2iso at a post-explosion time of 32. 1 d wit h 
the NTT observations on 14 Nov. 1999 of SN1999em faamuv et af]|200ll) . 
This observation corresp onds to an inferred post-explosion time of ~22 d 
iDessart & H illier 2006). These dates do not match. Instead, our choice is 
done so that the observations and models have the same color. Middle: Same 
as top, but now for the sl5el2iso at a post-explosion tim e of 75.8 d and the 
Keck observations on 5 Dec. 1999 ( Leonard et al.ll2002l) . corresponding to 
a post-explosion time of 42 d. Bottom: Same as top, but now for the model 
sl5el2iso at a post-explosion time of 306 d and the Lick observations of 26 
September 2000 (331 d after discovery). In each panel, the synthetic spectra 
have been scaled to correct for the magnitude mismatch shown in Fig. 1211 



5 Non-thermal processes will likely increase the electron density but this 
will only have a secondary effect on many of the lines since they are the 
coolants, and hence their strength is primarily determined by energy depo- 
sition from radioactive decay 
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are 1) the treatment of non-LTE for all level populations 2) the ex- 
plicit treatment of line-blanketing, 3) the consideration of all time- 
dependent terms appearing in the statistical-equilibrium, energy, 
and radiative-transfer equations, 4) the adoption of a hydrodynam- 
ical model describing the density/temperature/radius/velocity dis- 
tribution as well as the complex chemical stratification in the SN 
ejecta, and 5) the modelling of the entire ejecta at all times. The 
main approximations of our approach are spherical symmetry, ho- 
mologous expansion, and the treatment of radioactive decay as a 
local heating source (thus ignoring the effects of non-thermal exci- 
tation and ionisation). 

From our simulations we find/infer the following: 

(i) The photospheric evolution during the plateau phase shows 
a continuous radius increase up to 2.5-3 x 10 15 cm at ~80d; a de- 
crease and then a plateau in temperature at ~5000 K; and a contin- 
uous decrease in velocity as the ejecta optical thickness diminishes 
with expansion and recombination. This is in qualitative agreement 
with that inferred from observation, and is sim ilar to the trend we 
obtained for SN1987A jDessart & Hillierl2010l) . 

(ii) The synthetic bolometric luminosity shows a very gradual 
decline (not quite a plateau) during the photospheric phase. At 50 d 
after explosion our models have bolometric luminosities of ~ 6 x 
10 8 L0 for model sl5el2 and ~ 10 9 Lq for model s25el2. These 
are a factor of 2 to 3 larger than t he standard 50-d lumino sity of 
~3 x 10 s L for observed SNe II-P dBersten & Hamuy|2009h. Sim- 
ilar ov erestimates result from the simulations of lKasen & Wooslevl 
(2009), which ar e based on similar progenitor models. However 
the simulation o f lEastmanetaIUl994h . based on a more compact 
pre-SN RSG star, is in agreement with observations. This suggests 
that the pre-SN RSG stars produced by current stellar-evolutionary 
models have surface radii that are too large, on average. 

(iii) The density distribution of the photosphere is too steep in 
our models, at 10 to 30 d after core collapse, yielding line strengths 
that are weaker than observed prior to the onset of the recombi- 
nation phase. This may independently indicate that the progenitor 
surface radii are too large in our models. Steady state models in- 
dicate that density exponents of ~10 (i.e., p oc r -10 ) rather than 
~25 (valid at early times for the present models) are required to fit 
the line strengths. 

(iv) At the end of the plateau phase, our computed bolometric 
luminosities decrease by a factor of a few over a timescale of about 
a month until reaching a rate that corresponds to the instantaneous 
decay-energy production rate. In our simulations, the magnitude 
of the fading at the end of the photospheric phase is conditioned 
by two properties that are not obviously linked — the progenitor 
radius and the 56 Ni mass produced, and is thus expected to occur 
with a scatter in Nature. 

(v) Our multi-band light curves tend to show bell shaped curves, 
with a peak that occurs earlier for filters positioned at shorter wave- 
length, and a width that increases for filters positioned at longer 
wavelength. This morphology is similar to th at seen for SN1987A 
dPhillips et al1ll988l ; iDessart & HiUieril2010l) . but the larger RSG 
progenitors corresponding to our models causes the effective/color 
temperature to decline more slowly (and the more so for model 
s25e!2 because of it s larger-radius). Comparison with SN1999em 
dLeonard et afl [2002) confirms that decline in effective/color tem- 
perature is too slow, with the U B magnitudes remaining large for 
about 10-20 d. 

(vi) As for our simulations of SN1987A, we predict too strong 
UV and U -band fluxes during much of the photosphere phase. This 



is likely due to inadequate model atoms, or the H/He abundance 
ratio in the envelope, and is currently being investigated. 

(vii) In both simulations the photosphere resides in the outer 
1 Mq of the SN ejecta for ~ 50 d after the explosion. Thus spec- 
tra taken during this time do not reveal any information about the 
chemical and density structure of the bulk of the mass. The con- 
finement of the photosphere to the outer layers also provides a sim- 
ple explanation f or the low linear-polarisation measures of SNe II- 
P at early times (|Leon ard & Filippenko 2005 ; Leonar d~et alj|2006l : 
Chorn ock et al .120 id) , and the notion that early-time SNe II-P spec- 
tra are relatively unaffected by the details of the explosion (e.g., ge- 
ometry and 56 Ni production). Rather, SNe II-P evolve at constant 
composition through changes in ionisation alone — it is at such 
early-times that SNe II-P are best used for distance determinations 
based on spectroscopic modelling. Given the issues with the den- 
sity structure, and the confinement of the photosphere to the outer 
layers at early times, it is clear that higher resolution pre-SN models 
are needed to allow more accurate results in our present modelling 
approach - specifically we require that the outer SN ejecta be re- 
solved and accurately described all the way to large velocities. 

(viii) At the same color, photospheric spectra show excellent 
agreement with those of SN1999em. Nebular spectra, with the ex- 
ception of the H I spectrum, are also strikingly similar to those ob- 
served. The problem with the Hi lines is directly attributable to 
the neglect of excitation/ionisation by non-thermal electrons. In our 
models only forbidden lines of O, Na, Ca, Ti, and Fe, and resonance 
lines of Ca shine at a few 100 d after explosion (in the optical and/or 
near-IR ranges). At 1000 d after explosion, the spectra are domi- 
nated by the double line of Ol and numerous lines of Fel-II; no 
other species included in our calculations seems to contribute at 
such late times. 

(ix) The sl5el2 and s25el2 models show differences in peak 
luminosity and the rise-time to peak V-band brightness. However, 
at the same colour, photospheric spectra are very similar. This sug- 
gests that the composition contrast between the two progenitors is 
too small to be observable spectroscopically during the "plateau" 
phase. It also corroborates the spectral degeneracy of SNe II-P 
spectra and highlights the difficulty of distinguishing between dif- 
ferent progenitor masses or evolutionary paths based on "plateau" 
phase spectra. 

We find that nebular-phase spectra show differences that unam- 
biguously distinguish a lower-mass progenitor from a higher-mass 
progenitor in a standard SN II-P explosion. These are related to 
the property that stars with higher main-sequence masses have sys- 
tematically a higher helium-core mass, and shows up in late-time 
[O I] 6303-6363 A line profil es, which appear broader and stronger 
for a higher-mass progenitor (Dessart et al. 2Qlo3). 



In this work, we have demonstrated that there is much to learn 
from synthetic spectra of SNe II-P. Spectra contain additional infor- 
mation that is either absent or degenerate in multi-band light curves 
alone. Ultimately, it is through non-LTE time-dependent radiative- 
transfer simulations of a large set of non-rotating and rotating pre- 
SN RSG stars, evolved at various metallicities, that an accurate 
mapping between SNe II-P observations and progenitor/explosions 
properties will be possible. We plan to undertake this task in the 
near future. Securing a large dataset of photospheric- and nebular- 
phase observations of SNe II-P, covering all stages of temporal evo- 
lution, is essential for the fulfilment of this task. 
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